Mass loss and very low-metallicity stars 
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Abstract. Mass loss plays a dominant role in the evolution of massive stars at solar metallicity. 
After discussing different mass loss mechanisms and their metallicity dependence, we present the 
possibility of strong mass loss at very low metallicity. Our models at Z = 10~ 8 show that stars more 
massive than about 60 M may lose a significant fraction of their initial mass in the red supergiant 
phase. This mass loss is due to the surface enrichment in CNO elements via rotational and convective 
mixing. Our 85 M model ends its life as a fast rotating WO type Wolf-Rayet star. Therefore the 
models predict the existence of type Ic SNe and long and soft GRBs at very low metallicities. Such 
strong mass loss in the red supergiant phase or the iir-limit could prevent the most massive stars 
from ending as pair-creation supernovae. 

The very low metallicity models calculated are also very interesting from the nucleosynthesis 
point of view. Indeed, the wind of the massive star models can reproduce the CNO abundances 
of the most metal-poor carbon-rich star known to date, HE1327-2326. Finally, using chemical 
evolution models, we are able to reproduce the evolution of CNO elements as observed in the normal 
extremely metal poor stars. 
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INTRODUCTION 

The first generations of stars took part in the re-ionisation of the universe at the end of 
the dark ages (roughly 400 million years after the Big Bang). They are therefore tightly 
linked to the formation of the first structures in the universe and can provide valuable 
information on the early evolution of the universe. The first stars are thought to be more 
massive than solar metallicity stars dBromm & Larsonl . l2004l : [Schneider et all 2006|) and 



mass loss is expected to be very low at very low metallicities. The logical deduction 
from these two arguments is that a large fraction of the first stars were very massive at 
their death (> 100 M ) and therefore lead to the production of pair-creation supernovae 
(PCSNe). Unfortunately, the first mass ive stars died a long time ago and will probabl y 
never be detected directly (see however lScannapieco et alll2005l : lTornatore et alll2007|) . 



There are nevertheless indirect observational c onstraints on the first star s coming from 
observations of the most metal poor halo stars (|Beers & C hristlieb. 2005). These obser- 



vations do not sh ow the peculiar chemical signature of PCSNe (strong odd-even effects 



and low zinc, see He ger & WooslevL 12002). This probably means that at most only a few 
of these very massive stars (>100 M ) formed or that they lost a lot of mass even though 
their initial metal content was very low. In this paper, we discuss the possibility of strong 
mass loss at very low metallicities. 



MASS LOSS AND ITS DEPENDENCE ON METALLICITY 



At solar metallicity (Z ), mass loss has a crucial impact on the evolution of massive 
stars. It affects evolutionary tracks, lifetimes and surface abundances. It also determines 
the population of massive stars (number of stars in each Wolf-Rayet subtype for exam- 
ple). It influences the type of supernova at the death of the star (SNII, lb, Ic) and the final 
remnant (neutron star or blac k hole). Mass loss release s matter, containing newly pro- 
duced helium and metals (see IChieffi & LimongiL 2006, for the impact of mass loss on 
yields), and energy back into the interstellar medium in amounts comparable to super- 
novae (for stars above 30 M & ). Finally, it affects the hardness of the ionizing radiation 
coming from massive stars. It is therefore very important to understand mass loss in 
order to understand and model the evolution of massive stars. 

The metallicity (Z) dependence of mass loss rates is usually described using the 
formula: 

M(Z)=M(Z & )(Z/Z Q ) a (1) 

The expo nent a varies between 0.5-0.6 dKudritzki & PulsL l2000l : iKudritzkil l2002h and 
0.7-0.86 dVink et all EoOll: IVink & de KoterL l2005h for O-tvoe and WR stars respec- 



tively (See iMokiem et all . 20071 for a recent comparison between mass loss prescrip- 



tions and observed mass loss rates). Until very recently, most models use at best the 
total metal content present at the surface of the star to determine the mass loss rate. 
However, the surface chemical composition becomes very different from the solar mix- 
ture, due either to mass loss in the WR stage or by internal mixing (convection and 
rotation) after the main sequence. It is therefore important to kno w the contribution 
from each chemi c al spe cies to opacity and mass loss. Recent studies dVinket all 12000: 
Vink & de KoteiL |2005) show that iron is the dominant element concerning radiation 



line-driven mass loss for O-type and WR stars. In the case of WR stars, there is however 
a plateau at low metallicity due to the contributions from light elements like carbon, 
nitrogen and oxy gen (CNO). In the red supergiant (R SG) stage, the rates generally used 
are still those of iNieuwenhuijzen & de Jagerl dl990|) . More recent observations indicate 
that there is a very weak dependence of dust-driven mass loss on metallicity and that 
CNO elem ents and especia l ly nucleation seed component s like si licon and titanium are 



dominant dvan Loonl . l2000l 120061 : iFerrarotti & GailL l2006h . See Ivan Loon et ail d2QQ5h 



for recent mass loss rate prescriptions in the RSG stage. In particular, the ratio of carbon 
to oxygen is important to determine which kind of molecules and dusts form. If the ra- 
tio of carbon to oxygen is larger than one, then carbon-rich dust would form, and more 
likely drive a wind since th ey are more opaque than oxygen-rich dust at low metallicity 
dHofner & AndersenLl2007h . 

In between the hot and cool parts of the HR-diagram, mass loss is not well under- 
stood. Observations of the LBV stage indicate that several solar masses per year may 
be lost dSmith et alll2003|) and there is no indication of a metallicity dependence. Chro- 
mospheric activity could also pla y a role in stars haying s urface temperature similar to 
the Sun. Thermally driven winds (|Schroder & Cuntd . l2007|) and pulsations are still other 
ways to lose mass. 

The mass loss prescriptions used in the Geneva stellar evolution code are described 
in detail in iMeynet & Maederl d2005|) . In particular, the mass loss rates depend on 



metallicity as M ~ (Z/Z Q ) 0,5 , where Z is the mass fraction of heavy elements at the 
surface of the star. Models including the effects of both mass loss and especially rotation 
better reproduce the WR/O ratio and also the ratio of type Ib+Ic to type II supernova as 
a func tion of metallicity compared to non-rotating models, which underestimate these 
ratios (IMeynet & Maederu2005|) . Rotating models also better reproduce the ratio of blue 
to red supergia nts at low metallicities an d surface enrichments in nitrogen during the 
main sequence ( Maeder & Meynettl2000|) . 



FIRST STELLAR GENERATIONS 

As we saw in the previous sections, mass loss plays a crucial role in the evolution of 
solar metallicity stars. In this section, we discuss the importance of mass loss on the 
evolution of the first stellar generations. The first stellar generations are different from 
solar metallicity stars due to their low metal content or absence of it. First, very low-Z 
stars are more compact due to lower opacity. Second, metal free stars burn hydrogen 
in a core, which is denser and hotter due to the lack of initial CNO elements. This 
implies that the transition between core hydrogen and helium burning is much shorter 
and smoother. Furthermore, hydrogen burns via the pp-chain in shell burning. These 
differe nces make the metal fr ee (first) stars different from the second or later generation 
stars! (lEkstrom et all l2007h . Third, mass loss is metallicity dependent (at least for 



radiation-driven winds) and therefore is expected to become very weak at very low 
metallicity. Finally, the initial mass fun ction of the first s t ellar g enerations is expected to 
be top heavy below a certain threshold ( Bromm & LoebL 2003). 



Mass loss is expected to be very weak. What could change this expectation? In order 
to investigate this question, models at very low and zero metallicity were calculated 
with the same physics as in the models able to reproduce many observables around solar 
metallicity. These models were computed with an initial total angular momentum similar 
or slightly higher to the one contained in solar metallicity models with average initial 
rotational velocities of 300 km s . Models of metal free stars including the effect of 



rotation (jEkstrom et all 12006) show that stars may lose up to 10 % of their initial mass 
due to the star rotating at its critical limit (also called break-up limit). The mass loss 
due to the star reaching the critical limit is non-negligible but not important enough to 
change drastically the fate of the first generation stars. 



Th e situation is v ery different at very low but non-zero metallicity (jMeynet et al. . 



20061: iHirschiL 120071) . The total mass of an 85 M© model at Z = 10 8 is shown in Fig 



Q] by the top solid line. This model, like metal free models, loses around 5% of its 
initial mass when its surface reaches break-up velocities in the second part of the main 
sequence. At the end of core H-burning, the core contracts and the envelope expands, 
thus decreasing the surface velocity and its ratio to the critical velocity. The mass loss 
rate becomes very low again until the star crosses the HR diagram and reaches the RSG 
stage. At this point the convective envelope dredges up CNO elements to the surface 
increasing its overall metallicity. The total metallicity, Z, is used in this model (including 
CNO elements) for the metallicity dependence of the mass loss. Therefore depending on 
how much CNO is brought up to the surface, the mass loss becomes very large again. 
The CNO brought to the surface comes from primary C and O produced in He-burning. 



Rotational and convective mixing brings these elements into the H-burning shell. A large 
fraction of the C and O is then transformed into primary nitrogen via the CNO cycle. 
Additional convective and rotational mixing is necessary to bring the primary CNO to 
the surface of the star. 
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FIGURE 1. Structure evolution diagram of the 85 M Q model at Z = 10~ 8 . Coloured areas correspond to 
convective zones along the lagrangian mass coordinate as a function of the time left until the core collapse. 
The top solid line shows the total mass of the star. The burning stage abbreviations are given below the 
time axis. 

The strongest mass loss occurs in these models in the cooler part of the HR diagram. 
Dust-driven winds appear to be metallicity independent as long as C-rich dust can form. 
For this to occur, the surface effective temperature needs to be low enough (log(T e ff) < 
3.6) and carbon needs to be more abundant than oxygen. Note that nucleation seeds 
(probably involving titanium) are still necessary to form C-rich dust. It is not clear 
whether extremely low-Z stars will reach such low effective temperatures. This depends 
on the opacity and the opacity tables used in our calculations did not a ccount for the 
non-standard mixture of metals (high CNO and low iron abundance, see lMarigoLl2002l 
for possible effects). It is interesting to note that the wind of the 85 M Q model is richer 
in carbon than oxygen, thus allowing C-rich dust to form if nucleation seeds are present. 
There may also be other important types of wind, like Chromospheric activity-driven, 
pulsation-driven, thermally-driven or continuum-driven winds. 



NUCLEOSYNTHESIS AND CHEMICAL EVOLUTION 
The most metal poor star known to date, HE 1327-2326 



Significant mass loss in very low-Z massive stars offers an interesting explanation for 
the strong enrichment in CNO elements of the most met al poor stars observed in the 
halo of the galaxy (see iMeynet et al. . 20061: iHirschiL 2007 ). The most metal poor stars 
known to date, HE 1327-2326 (|Frebel et all l2006h is characterised by very high N, C 
and O abundances, high Na, Mg and Al abundances, a weak s-process enrichment and 
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FIGURE 2. Composition in [X/Fe] of the stellar wind for the Z = 10 models (solid lines). For 
HE1327-2326 (red stars), the best fit for the CN O elements is obtained by diluting the composition of 
the wind of the 40 M model by a factor 600 (see iHirschil 12007, for more details). 



depleted lithium. The star is not evolved so has not had time to bring self-produced 
CNO elements to its surface a nd is most likely a sub gia nt. By using one or a few SNe 
and using a very large mass cut. lLimongi et al.l (|2003h and llwamoto et al. (2005) are able 
to reproduce the abundance of most elements. However they are not able to reproduce 
the nitrogen surface abundance of HE 1327-2326 without rotational mixing. A lot of the 
features of this star are similar to the properties of the stellar winds of very metal poor 
rotating stars. HE1327-2326 may therefore have formed from gas, which was mainly 
enriched by stellar winds of rotating very low metallicity stars. In this scenario, a first 
generation of stars (PopIII) pollu tes the interste llar medium to very low metallicities 
([Fe/H]— 6). Then a PopII.5 star dHirschil . l2005h like the 40 M & model calculated here 
pollutes (mainly through its wind) the interstellar medium out of which HE1327-2326 
forms. This would mean that HE 1327-2326 is a third generation star. In this scenario, the 
CNO abundances are well reproduce d, in particular th at of nitrogen, which according to 
the new values for a subgiant from lFrebel et alJ (|2006h is 0.9 dex higher in [X/Fe] than 
oxygen. This is shown in Fig. [2] where the abundances of HE1327-2326 are represented 
by the red stars and the best fit is obtained by diluting the composition of the wind of 
the 40 Mq model by a factor 600. When the SN contribution is added, the [X/Fe] ratio 
is usually lower for nitrogen than for oxygen. It is interesting to note that the very high 
CNO yields of the 40 M Q stars brings the total me tallicity Z above the limit for low mass 
star formation obtained in lBromm & Loebl (|2003r) . 



Primary nitrogen 



One of the most stringent observational constraint at very low Z is a very high primary 
N production. This requires extremely high primary 14 N production in massive stars, 
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FIGURE 3. Chemical evolution model predictions of the N/O and C/O evolution, in the galactic halo, 
for different stellar evolution inputs. The solid curves show the predictions of a model without fast rotators 
at low metallicities. The dashed and dotted lines, almost overlapping, show the effect of including a 
populati on of fast rotato r s at lo w metallicities (the dotted line includes also the Z=0 fast rotators). For the 
data see Chiapp ini et al.l J2006) and references therein. 



about 0.1 M Q per star (~0.15 M© used in the heuristic model offchiappini et al] . l20"05j) . 
Upon the inclusion of the new stellar calculations of iHirschil (|2007l) for Z=10 8 in 
a chemical evolution model for the galactic halo with infall and outflow, both high 
N/O and C/O ratios are obtained in the very-metal poor metallicity range in agreement 
with observations (see details in Chiappini et all l2006b|) . This model is shown in Fig. 
[3] (dashed magenta curve). In the same figure, a model computed without fast rotators 
(solid black curve) is also shown. Fast rotation enhances the nitrogen production by ~3 
orders of magnitude. These results also offer a natural explanation for the large scatter 
observed in the N/O abundance ratio of normal metal-poor halo stars: given the strong 
dependency of the nitrogen yields on the rotational velocity of the star, we expect a 
scatter in the N/O ratio which could be the consequence of the distribution of the stellar 
rotational velocities as a function of metallicity. Furthermore, the strong production of 
primary nitrogen is linked to a very active H-burning shell which contributes a large 
part of the total luminosity of the star. Hence, the energy produced by the helium core 
is reduced, reducing the efficiency of the 12 C(a,y) 16 reaction. As a consequence, 
less carbon is turned into oxygen, producing high C/O ratios. Although the abundance 
data for C/O is s till very uncerta in, a C/O upturn at low metallicities is suggested by 
observations (see lAsplund , 2005, and references therein). 



GAMMA-RAY BURSTS AND PAIR-CREATION SUPERNOVAE 



Long and soft gamma-ray bursts (GRBs) have now been firmly connected to the death 
of type Ic supernovae (see lWoosley & Blooml 20061 fo r a rec ent review). In one of the 
most promising models, the collapsar model ( Woosleyl Il993|) . GRB progenitors must 



form a black hole, lose their hydrogen rich envelope (become a WR) and retain enough 
angular momentum in their core during the pre-supernova stages. The strong mass loss 
discussed in the previous section makes it possible for single massive stars in the first 
stellar generation s to become W R stars and even to retain enough angular momentum to 
produce a GRB ( Hirschil 2007 ). A wide grid of models around Z Q shows that the fast 
rotating WO stars are produced only at metallicities equ al to or lower than that of the 
LMC, in agreement with observations (IStanek et all 12006). These models predict that, at 
low Z, one GRB occurs for every 10 core co llapse supernovae, w hich is in the upper limit 



of what is allowed by the observations (see lHirschi et all I2005L for more det ails). More 



recent models including the effects of magnetic fields (|Yoon & Langeii l2005|) show that 
another mechanism is possible to produce GRBs at low Z. This mechanism is the quasi- 
homogeneous evolution of very fast rotating massive stars. In this scenario, WR stars are 
produced by mixing and not mass loss. This last scenario however does not predict GRB 
at metallicities equal or higher than th e SMC. This upper limit is too low compared to 
recent observations (IStanek et all 120061). Taking into account the anisotropy in the wind 
of these fast rotating stars (IMeynet & Maedeii 120071) may help reduce the discrepancy 
between models and observations. Note that the downward revision of solar metallicity 
( Asplundi 2005) may help resolve the problem. 

Apart from GRBs, pair-creation supernovae (PCS Ne) are very energetic ex plosions, 
which could be observed up to very high redshifts (IScannapieco et al. . l2005h . PCSNe 
are expected to follow the death of stars in the mass range between 100 and 250 Mq, 
assuming that they do not lose a significant fraction of their mass during the pre- 
supernova stages. Amongst the very first stars formed in the Universe, one expects 
to have PCSN due to the lack of mass loss and to the low opacity unable to stop the 
accretion on the star during its formation. However, the EMP stars observed in the halo 
of the gal axy do not show the pecu liar chemical signature of PCSN (strong odd-even 
effect, see iHeger & WooslevL [2002). This means that either too few or even no PCSN 
existed. One possible explanation to avoid the production of very low-Z or metal free 
PCSNe is the strong mass loss in the cool part of the HR diagram due to the surface 
enrichment in CNO elements induced by r otational and convec tive mixing (see previous 
section) or the star reaching the ^r-limit (lEkstrom et alll2007|) . 



CONCLUSION 

For the most massive models (M > 60 Mq), significant mass loss occurs during the red 
supergiant stage assuming that CNO elements are important contributors to mass loss. 
This mass loss is due to the surface enrichment in CNO elements via rotational and 
convective mixing. The models predict the production of WR stars for an initial mass 
higher than 60 M Q at Z = 10~ 8 and our 85 M & model with V[ n i = 800km s~ ! becomes 
a WO. Therefore SNe of type lb and Ic are predicted from single massive stars at these 
low metallicities. The 85 Mq model retains enough angular momentum to produce a 
GRB but the calculations did not include the effects of magnetic fields. 

The stellar yields were calculated for light elements. These yields were used in a 
galactic chemical evolution model and successfully reproduce the early evolution of 
CNO elements ( Chiappini et all 120061) . Finally, a scenario is proposed to explain the 



CNO abundances of the most metal poor star known to date HE1327-2326. 
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